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A generalized autocatalytic model for chiral polymerization is investigated in detail. Apart from 
enantiomeric cross-inhibition, the model allows for the autogenic (non-catalytic) formation of left 
and right-handed monomers from a substrate with reaction rates el and er, respectively. The 
spatiotemporal evolution of the net chiral asymmetry is studied for models with several values of 
the maximum polymer length, N. For N — 2, we study the validity of the adiabatic approximation 
often cited in the literature. We show that the approximation obtains the correct equilibrium values 
of the net chirality, but fails to reproduce the short time behavior. We show also that the autogenic 
term in the full N = 2 model behaves as a control parameter in a chiral symmetry-breaking phase 
transition leading to full homochirality from racemic initial conditions. We study the dynamics of 
the N ^ oo model with symmetric {el ~ er) autogenic formation, showing that it only achieves 
homochirality for e < Ec, where Ec is an A''-dependent critical value. For e < Ec we investigate 
the behavior of models with several values of A'^, showing that the net chiral asymmetry grows as 
tanh(A). We show that for a given symmetric autogenic reaction rate, the net chirality and the 
concentrations of chirally pure polymers increase with the maximum polymer length in the model. 
We briefly discuss the consequences of our results for the development of homochirality in prebiotic 
Earth and possible experimental verification of our findings. 



I. INTRODUCTION 

Life is dependent upon biomolecular spatial asymmetry. Even though the origin of this asymmetry remains unknown 
[7], it is frequently argued that hfe could neither exist nor originate without molecular chirahty [2,. Under this view, 
the question of the origin of homochirahty is intrinsic to that of the origin of life: if we want to understand how 
life first emerged on Earth (or possibly elsewhere), we must understand how homochirality emerges from basic chiral 
building blocks (enantiomers) . 

The mystery here may not be so much which chirality but why chirality: although the sugars found in DNA and RNA 
are right-handed and protein amino acids are left-handed, laboratory syntheses yield racemic mixtures A pertinent 
exception is the reaction studied by Soai and coworkers, which demonstrates bifurcation to a single handedness |23j . 
Although the Soai reaction has probably little relevance to life on early Earth, it exhibits the features required of a 
successful model: autocatalysis and enantiomeric cross-inhibition fed by enzymatic enhancement. 

Inspired by Frank's pioneering work [H], Sandars recently proposed an autocatalytic model including enantiomeric 
cross-inhibition [2T]. Sandars' model has been further investigated in work by Brandenburg and collaborators OH], 
Wattis and Coveney |21I, and Gleiser and Thorarinson [9l|T0]. One of the model's shortcomings is that enantiomeric 
feedback is generated by the longest polymer chain in the model. This requires somewhat unrealistic initial conditions, 
as they must include initial asymmetric concentrations of the maximum length polymers with no intermediate length 
polymers present in the reactor pool (i.e., requiring the presence of Ljv and Rn with some asymmetry and not 
requiring the presence of L„ or i?„ for n < N). In a more realistic bottom- up approach, initial conditions should 
depend only on the substrate and possibly on a very small asymmetry between the two types of monomers, perhaps 
caused by random fluctuations in concentration [S]. At an even more basic level, one might desire the system to start 
with no polymers, monomers, and even no substantial concentrations of substrate, so long as autogenic production of 
monomers is active. Such non-catalytic production terms for monomers of both chiralities were briefly discussed by 
Wattis and Coveney Saito and Hyuga [20], and by Brandenburg et al. [2111] but have not yet been thoroughly 
explored. 

In the present work, we extend Sandars' model by allowing for spontaneous monomer production from the substrate 
and for spatial dependence of the concentrations. We study in detail the spatiotemporal dynamics of the reaction 
network for several values of maximum polymer length N, including the iV — > oo limit. For the N = 2 (or reduced) 
model, we investigate the validity of the adiabatic approximations employed by Brandenburg and Multamaki 0], 
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Gleiser and Thorarinson [5], and by Gleiser [TO]. We show that the introduction of autogenic monomer production 
influences the evolution of the net chiraUty by inducing a chiral symmetry-breaking phase transition. For models with 
> 2, we obtain the net chirahty as a function of N, showing that it is sensitive to the autogenic reaction rates, 
^L(i?)i for left (right )-handed monomers. For symmetric autogenic reactions, (eL(fl) = e), we show that above a critical 
value £c no net chirality is produced: the racemizing autogenic term overwhelms the chirality-inducing enantiomeric 
cross-inhibition . 

The paper is organized as follows. In sectionjlT] The General Model, we introduce the polymerization model and the 
relevant rate equations. The key features of autocatalysis, enantiomeric cross-inhibition, and enzymatic enhancement 



are discussed. In section III The Reduced Model, we describe the reduced {N = 2) model and the interpretation of 
the net chirality as an order parameter satisfying an effective potential. We also test the validity of the adiabatic 
approximations and describe the dynamics of a chiral symmetry-breaking transition. In section |IV[ The N —^ oo 
Model, we investigate the limit iV — > oo. For symmetric autogenic reactions, we obtain the phase diagram for the net 
chirality as a function of e, showing that above a critical value Sc homochirality is unattainable. We also compute 
the average polymer length as a function of e, showing that as £ ^ Ec, large chiral chains are strongly suppressed. 
We then examine the spatiotemporal dynamics for several models with iV < 10 and iV ^ cxd in section [V] Numerical 
Results: Maximum Polymer Length and Chiral Asymmetry. We obtain an expression for the net chirality as a function 
of maximum polymer length iV as £ varies. We also demonstrate that larger polymers achieve higher chiral purity. In 
section [VTl Conclusion, we present a brief summary of our results. 



II. THE GENERAL MODEL 



Sandars [2T] proposed a model for prebiotic homochirality whereby homochiral polymers evolve from a gradual build 
up of chiral monomers. The model includes the key features for chiral ampliflcation, enantiomeric cross-inhibition 
and enzymatic enhancement. In this section, we will present our generalization of Sandar's model, and discuss how it 
compares with other models proposed recently in the literature. Given that little is known of the specifics of prebiotic 
chemistry, it is important to examine the effects of plausible features expected to be present in realistic models. 



A. Modeling Polymerization 

Consider a chirally pure chain of n left-handed monomers, L„. It can grow on either end by addition of a left 
or right-handed monomer, Li or Ri, with rate coefficients fcs £^nd kj, respectively. The attachment of the "wrong" 
monomer, in this case a right-handed monomer, inhibits that end of the polymer from further growth. This process is 
referred to as enantiomeric cross-inhibition: addition of a monomer of the opposite chirality to one end of a growing 
chain terminates further growth on that end of the chain |!l5, . The reaction network may be written as 

Ln + ii — > Ln+l, 

Ln + Rl ^ LnRl, 
Ll + LnRl Ln^iRi, 
Rl + LnRl RlLnRl, 

(1) 

supplemented by the complementary reactions obtained by exchanging L ^ R. 

In addition, the reaction network includes a substrate S from which monomers of both chiralities are generated. The 
substrate may generate monomers spontaneously or through enzymatic enhancement by already existing polymers. 
The rate equations governing these processes are 

O — > L/i, J — > L/i, o — > -Ttl, 

S'JiRi, s'^^'^Ri, S'^^'-Li, (2) 

where El and er are the rate coefficients for direct (autogenic or non-catalytic) production of monomers from the sub- 
strate, kc is the rate coefficient for enantiomeric feedback, Cl{Cj^) denote the enzymatic enhancement of left(right)- 
handed monomers, and the coefficients p and q are given by, p — ^(1 + /) and q ~ ^(1 — /), where / is the fidelity of 
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the enzymatic reactions. Note that for complete fidehty (/ = 1), g = 0, and the two rightmost processes in eq. [2] do 
not occur. 

For simplicity, Sandars chose = and Cr = R^, where N is the maximum polymer length. However, given 
that the exact functional dependence is not known, alternative suggestions have been made, such as Cl ~ "^^n 
and Cr — 'Y^nRn [3J [23] (henceforth WC and BAHN, respectively). Here, we allow for polymers of all lengths to 
participate equally in the autocatalytic feedback and choose Ch — ^ Ln and Cr = "^Rn- In contrast to the case 
where the enzymatic enhancement is dependent only on the longest iV-length polymers, this choice allows homochiral 
polymers of all lengths to participate in the autocatalytic feedback and removes any dependence on the specific cutoff 
value of N. This choice also allows the network to start with initial conditions which do not require the artificial 
presence of small, left-right asymmetric concentrations of the highest order polymer, a limitation of Sandar's model. 
In a bottom-up approach, long chains should be products of a gradual build-up: it is unlikely that they would be 
present ab initio. 

The rate coefficients El and sr describe the production of monomers from substrate without the need for catalytic 
feedback. Given that little is known of the details of prebiotic chemistry, it is important to consider the effects of such 
terms on the reaction kinetics. Furthermore, chiral monomers could also have been fed to the early Earth environment 
from outer-space, as supported from findings of chiral amino acids in, e.g., the Murchinson meteorite ^ l19|. Similar 
terms have previously been introduced by WC, Saito and Hyuga [20] (henceforth SH), and by BAHN, although WC 
and SH did not allow for biased production. In the model of WC, a single non-catalytic production term e = Sl = £r 
is introduced. WC consider only e <C 1 to allow for simplified initial conditions (i.e., starting with no substrate, 
monomers, or polymers present but with a continuous source term for the substrate) and did not study the effect of 
this term on the dynamics of the model or its equilibrium solutions, as we do here. BAHN introduced the non-catalytic 
production terms Cql and Cqr. Comparing to our model, one finds — kcCo^ and Sr = kcC^R. BAHN found that 
the introduction of finite Cql or Cqr with Cql ^ Cqr leads to imperfect but still efficient bifurcation in the feedback 
fidelity and to the absence of a racemic equilibrium solution. These authors suggest that homochirality is primarily 
the result of the instability of the racemic solutions which is barely modified by the presence of finite Cql or Cqr. As 
a consequence, most of BAHN's work assumed Cql — Cqr — 0. Here, we find that even in the case = Er — e, the 
presence of these terms can, in fact, infiuence the instability of the racemic solutions: for e > ec{N), homochirality 
is impossible (as will be discussed, Sc{N) is a critical value which varies with maximum polymer length). A similar 
conclusion was reached by SH, who studied autocatalytic reactions in a simplified system with two chiral species. 
SH found that random production of chiral molecules from non-catalytic processes must occur to generate initial 
concentrations of chiral products and that there exists a critical value for the non-catalytic production rate: above the 
critical value, the system can obtain only fractional homochirality. We arrive at a qualitatively stronger conclusion: 
ec{N) varies with the maximum length of polymers allowed in the model, N, and the net chirality for a given e < Sc 
increases with increasing N. 

Gleiser \10i investigated the case where the bias originated in the enantiomeric enhancement terms through the 
definitions Ql = kc{l + g/2)[S]CL and Qr — fcc(l — 9/'^)[S]Cr, where g <C 1. Exploring the efficacy of both 
intrinsic and extrinsic infiuences through the parameter g, it was found that in the diffusive regime only a fairly large 
bias (g > 10"^) could influence the evolution of reasonably-sized reactor pools toward homochirality in early Earth. 
In this regime, small sources of chiral symmetry breaking, such as parity violation in the weak interactions where 
g ^ 10-i^[16ll25], were ruled out. 

With the above processes of deriving monomers (both non-catalytically and through enzymatic enhancement) 
present in the model, it is possible to start from racemic initial conditions and, through chiral amplification, arrive at 
a homochiral final state. One can imagine initially having zero concentrations of polymers, monomers, and substrate 
with the substrate being added continuously to the system starting at t = 0. As pointed out by WC, if £L(fl) = no 
monomers are produced and polymerization is never initiated. (Models without autogenic production usually start 
with an unrealistic nonzero concentration of large polymers.) We therefore require the presence of non-zero autogenic 
production terms, el and £r, to produce the initial concentrations of monomers. In essence, they jump-start the 
reaction network. 

There are two ways in which homochirality may be achieved once there is a nonzero concentration of substrate and 
possibly minute initial concentrations of monomers. If an explicit bias is introduced such that El ^ £r, the system 
requires an initial concentration of only the substrate in order for chiral amplification toward the favored chirality 
to occur. If £l = er, a small bias must be introduced between the initial concentrations of left and right-handed 
monomers. The needed asymmetry may be the result of random fluctuations in concentrations of large numbers of 
molecules [5], or it may be due to environmental disturbances [SHU]. We will explore both these cases, El = er and 
£l 7^ in the context of the present paper. 
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B. Polymerization Equations 

Taking the above into consideration, and accounting for the losses of Li and Ri to lengthen growing chains, the set 
of reaction rate equations governing the various concentrations of a A''-polymer network, are 

"^'^"^ = 2fcs[ii][in-i]-2[L„](fcs[ii] + fc/[i?i]), (3) 



dt 

d[R, 



dt 

d[L„Ri] 
dt 



= 2ks[Rl][Rn-l] - 2[Rr,]{ks[Rl] + ki[L,]), 

= ks[Li][Lr,-iRi] + 2fc/[i?i][L„] - [L„Ri]{ks[Li] + ki[Ri]), 



^t^^ = ks[Ri][Rn-iLi]+2ki[L,][Rn]-[RnL,]{ks[Ri] + ki[Li]), 

where 2 < n < N. We note that WC imposed the condition that [-£/„-Ri] and [-Rnii] were inhibited from further 

growth by enantiomeric inhibition, so the last two eqs. above were not included in their model. Here, we assume that 
these chains can still grow on the opposite end from that attached to a "wrong" enantiomer. For n = 2, we have the 
equations, 

= 2ki[Ry][L^] - [L2Ri]{ks[L^] + ki[R^]), (4) 
= 2ki[L^][R^]-[R^L^]{ks[Ri\ + ki[L^]). 

at 

In the above equations, factors of 2 arise when a monomer can attach to either end of a growing polymer. For 
dimers, one must discount this factor to account for the interaction of two single monomers. Thus, we must write 

d[L2] _ , rr 12 



dt 
d[R2] 



dt 

The evolution equations for left and right-handed monomers are, 

d[L 



= ks\LiY -2\L2]{ks[Lr] + ki[Rx]), (5) 
= ks[Rif - 2[R2]{ks[Ri] + 



dt 
d[R 



- = sl[S] + Ql-Xl[Li], (6) 



^^ii = eR[S]+QR-XR[Ri], 
where 

7V-1 N N-1 N 

Al = 2ks [Lr,] + 2ki + fcs ^ [LnR\ + ki 

n=l n=l n=2 n=2 

N-l N N-l N 

Xr = 2ks [Rn] + 2ki Yi^n] +ksYl i^r^L] + kj ^[L„i?]. 
n=l n=l n=2 n=2 

Xl and Xr quantify the losses from monomer populations associated with attachment of a monomer to a growing 
polymer, el and er are rate coefficients for the direct production of monomers from the substrate. The choice of 
N is model dependent. In the model of Sandars, N = N — 1, and in the model of BAHN, N = N. For the sake of 
simplicity, and to be consistent with other terms in the decay rates Xl and A^, we choose N = N — 1. 
Assuming that the substrate is maintained by a source, Qs, its evolution equation takes the form 

^ = Qs~{eL + eR.)[S]-{QL + QR) (7) 
where Ql and Qr are the source terms for production of monomers from enzymatic enhancement and take the form 

QL = kc[S]{pCL + qCR), (8) 

QR = kc[S]{pCR + qCL), (9) 
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where p and q are r elated to the fidehty as defined in section II A In this paper we set / = 1 (i.e., q = and p = 1). 
As stated in section 



II A 



we choose Cl = Y.n=2\-^ri] and Cr = Y.^^2{^n] 



This completes the model. However, note that the rate equation for [Lii?i], which may be found in the polymer- 
ization model of Sandars, has been left out. Due to the effects of enantiomeric cross-inhibition, the further growth of 
this dimer is terminated; consequently, it does not affect the essential dynamics of the model and we have chosen not 
to include it since it has no bearing on the net chirality. Even so, it is important to acknowledge that hybrid dimers 
of this form are being created and thus are present in the system. One would need some method of removing these 
dimers, as well as longer chains containing the minority chirality from the reactor pool. Sandars suggested a recycling 
process involving a non-chiral intermediary to remove initial quantities of the unfavored enantiomer |22j . We also 
note that even though most present-day biomolecules polymerize only unidirectionally fl8', given that we don't know 
the conditions of prebiotic chemistry, we chose to consider the more generic case where monomers can attach to both 
ends of a growing polymer (see eqjs]). The ensuing kinetics should be qualitative similar in the two cases, as the main 
differences will be that certain rate coefficients will be halved and that the equations for [L„i?i] and [i?„Li] are no 
longer necessary as these polymers do not grow. 

In order to discuss the net chirality in later sections, we introduce the following variables of the enantiomeric excess 

§ = ~ g = J2Ln-J2 ^ ^ Ln - Rn ^ ^^^y^ 

Ll + Rl J2-^n + J2 Ln + Rn ~ 

which describe different measures of chiral purity of homochiral polymers. In particular, 9 accounts for the net chiral 
excess from all homochiral chains. A special case occurs in the truncated N = 2 model with adiabatic approximations 
(discussed in the next section) , where 5 is the only measure we may use. The parameter ry„ describes the net chirality 
for chains of a specific length n. In the case n = N , this is equivalent to the measure of enantiomeric excess used 
by Sandars |21J. In the following sections, we will compare the temporal and spatiotemporal evolution of systems for 
various values of A^: N — 2, 5, 10, and in the limit N ^ oo. We start by considering the simplest case, the N ^ 2, 
or reduced, model. 



III. THE REDUCED MODEL 



We first consider the truncated N — 2 model of Brandenburg and Multamaki |3] (henceforth BM), where the longest 
polymers formed are dimers. Apart from its simplicity, the attractiveness of this model relies on the fact that dimers 
play a similar catalytic role to that demonstrated in the Soai reaction, We will examine if such a truncation 
maintains the essential dynamics of a higher N model. 



A. The Reduced Model with Adiabatic Approximation 



BM simplified the reaction network for the N — 2 system further by assuming that the rate of change of the 
concentrations of the substrate, [S], and of dimers, [L2] and [-R2], arc much slower than that of the monomers, [Li] 
and [Rl]. This approximation is known as the adiabatic elimination of rapidly adjusting variables [14, . We start by 
testing the validity of this approximation. 

It proves convenient to introduce dimensionless symmetric and asymmetric variables, S = X + Y and A = X — Y , 
respectively, where X = [Li]{2ks /QsV^^ and Y = [Ri]{2ks/Qsy^^- For ks/ki = 1, after some algebra, the 
polymerization equations reduce to 

A,-'f ^ l-5». (U) 

idA ^ 2S{-fL-jR + A) 
''o dt 2(7^-1-7^)5 + 52 + ^2 

where the parameter Aq = {2ksQs)^^^ has dimensions of inverse time, and we have introduced the parameters 

7L = 2( — )^— , 7i? = 2( — )2— . (12) 
Qs kc Qs kc 



For 'jR = jL = 0, the above equations reduce to the models of BM and of Gleiser and Thorarinson f3]. 
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FIG. 1: Left: Potential V{A) for = = "/ for several values of 7. Right; Potential for 7_r 7^ 7_l with varying asymmetry 
between the two parameters. In both graphs, the thick line corresponds to 7_r = 7_t = 7 = 0. 



The fixed point of the first equation in eqs. [TT]is easily found to be 5 = 1: the system tends quickly toward this 



value at time-scales of order Aq . With 5=1, the equation for the chiral asymmetry becomes 

^_^dA _ A-A^ + 2^L{l-A)-2-fR{l + A) 



dt 



Writing the above in the more suggestive form, 
potential 



A^ 



(13) 



-dV/dA (the dot denotes time derivative), one obtains the 



V{A) = ^ - ln[l +A^ + 2(7i + 7^)] 



2(7L - 1r) 
v/1 + 2(7i + 7;^) 



arctan 



A 



Vl + 2(7L+7j^) 



(14) 



The potential V{A) is symmetric under A —A for 7^ = 7^^ = 7, and asymmetric for 7^ 7^ 7/j. Figure [T] shows 
y(^) for the symmetric and asymmetric cases. 

Considering the symmetric case, 7 = 7l = Jr, the potential has fixed points a,t A — and A = ±-y/l — 47. The 
point ^ = describes the racemic solution and is unstable when the chiral fixed points at ^ = ±^1 — 47 are real. 
The enantiomeric excess is given by 



A 



±v/l-47. 



(15) 



Clearly, an enantiomeric excess is possible only for 7 < 0.25, and complete homochirality (5 = ±1) occurs only for 
7 = 0. For 7 > 7c = 0.25 a racemate {A = 0) is the only stable solution. We see that 7 plays a similar quantitative 
role as the temperature in systems in the Ising universality class [12J. We will further explore this analogy below. 

As discussed previously, BAHN found that in cases with an asymmetry between non-catalytic production terms, 
there is no longer a perfectly racemic stable solution. For 7l(7_r) large enough, the left(right)-handed chirality 
invariably will achieve dominance (c/. figure [l] right). For the N = 2 truncated model, the asymmetry between 7/, 
and 7fl' plays a similar qualitative role to the biasing parameter g introduced by Gleiser [10|. 



B. Chiral Selection as a Phase Transition 



The form of the potential in eq. [14] suggests the possibility of chiral symmetry breaking initiated by a temporal 
or spatiotemporal varying parameter ^{t) or 7(x, t), respectively. We treat spatial evolution in the following sec- 
tions. Here, we investigate the parameter 7 and its qualitatively similar role to that played by temperature in phase 
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transitions (see [H]). We consider the symmetric case ^ = = 1r and rewrite eqs. 12 as 



7-2(|^)^f, (16) 
Qs kc 

where e = El = £r- 

Consider a temporally varying parameter 7(t). One can imagine this to occur with a change in temperature, 
pressure, or concentration of reactants as the system evolves. Starting with a large 7, for instance 'y{t = 0) = 1, the 
production rate of monomers of both chiralities is relatively high (i.e., s > kc, holding ks fixed). As 7 decreases, the 
rate of direct production of monomers from the substrate slows, and enantiomeric feedback, governed by the reaction 
rate kc, becomes the dominant mechanism for formation of monomers (i.e., when s < kc)- As studied extensively in 
previous work [3J 1211 124] . enantiomeric feedback is the driving force that leads to eventual chiral bifurcation. Thus 
the critical point in this model occurs when the dominant mechanism for monomer formation switches from autogenic 
production from the substrate to enantiomeric feedback. The transition occurs at 7 = 7c (7c = 0.25 for the symmetric 
potential in figure [T| and complete chiral separation occurs as 7 ^ 0. 

In this scenario, as "f(t) passes through 7^, the system undergoes a phase transition from a racemic state to a chiral 
state. One can think in analogy with a ferromagnet: as it is cooled through the Curie point, it settles into one of the 
two degenerate ground states (magnetization directions) . The choice is determined by fluctuations about equilibrium. 
In the absence of an external magnetic field (the biasing factor, here equivalent to choosing 7^ ^ 7^), the choice is 
random; there is an equal probability that either ground state will be chosen. The net chirality plays the role of the 
net magnetization and we observe domains of both chiralities emerge as 7 is lowered through 7^. 

If one then considers the parameter 7 to be both spatially and temporally dependent, 7 = 7(2;, t) - a realistic 
possibility given that, in general, the temperature, pressure, or concentrations are likely to be inhomogeneous in early 
planetary environments - the phase transition can occur at different places at different times. For example, if 7 is 
dependent on temperature, then j(T) will be nonuniform, as the temperature itself may be a function of space and 
time (r(x, t)). In this case, the phase transition may not occur everywhere at once and chirality will not be uniform 
throughout space. Chiral domains will then compete for dominance as discussed in Gleiser |10| . The formation of 
chiral domains for a model phase transition will be illustrated below in the context of the full N = 2 model. 



C. The Reduced Model Without Adiabatic Approximation 

To lay out the groundwork for larger N, we explicitly write out the equations for the full N — 2 model without 
adiabatic approximations. This will also allow us to test their validity. We define the dimcnsionless symmetric and 
asymmetric variables, 

51 = {^)H[Li] + [i?i]), A, = i^)H[L^] - [i?i]), 

Qs Qs 

52 = {^)H[L2] + [R2]), A2 = (^)^([i2] - [R2]), 

Qs Qs 

describing the evolution of the concentrations of monomers and dimers, and the dimensionless variable ip — {^f)^ [S] 
describing the evolution of the substrate. The network equations are then (setting kj/ks = 1), 

Aq^*^ = l-^'i{lL + lR)lp- nijS2, 

Ao'^ = ^(^\{il + 1r)+S2^^-SI 

Ao'^ = ^i(^l{lL~lR)+A2^ij-SiAi, 

Ao^^ - \iS! + Al)-S2S,, 



Ao'^ = ^5i^i-^25i, (17) 



where, as before, Aq — {2ksQsY^'^ ■ The parameters 7^ and 77? were defined in eq. 
parameter k = 



12 



and we introduced the new 
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FIG. 2: Time evolution of the net chirality for N = 2 models with (continuous line) and without (dashed line) adiabatic 
approximations. Note that the asymptotic values of S are equal, and given by eq. |15[ 

For the remainder of this work we assume kj = ks = kc {k = 0.5): the rates of attachment of a monomer of the 
same or opposite chirality are equal, and these rates are in turn equal to the rate of catalytic production of monomers 
from the substrate by enantiomeric feedback. The other relevant rates that appear in our reaction network, el 
and Eft, which describe non-catalytic monomer production, should be much slower [31 [23]. For example, WC chose 
£ w 10~^s~^. To show that the condition e ^ 1 is satisfied, we use the nominal values of ks ~ lO^^^cm'^s"^ and 
Qs ~ lO^^cm^^s^^ with kc = ks to find 

7i(fl,) ^2V2x lO^ELiR) (18) 
Noting that 7 should be dimensionless, we obtain eL(fl) ^ lO^^s^^ for < 7 < 

D. Comparison of Reduced Models: Testing the Adiabatic Approximation 

In the previous subsections we introduced the N = 2 model with and without adiabatic approximations. Since the 
adiabatic approximations have been used in several models [4, |9j [10] , it is important to test their validity. In figure 
[2] we show the temporal evolution of 6 for both models with 7 = 0.1. The initial conditions were the same in both 
models: an asymmetry ([Li] — [i?i])|t=o = 10~^ in monomer concentrations. All ordinary differential equations in this 
work were solved numerically using Mathematica. 

If, as suggested by BM, adiabatic elimination does not affect the accuracy of steady-state solutions (where dS/dt 
as i — > 00), we should expect the asymptotic values to be equivalent. In figure |2] we see that the asymptotic solutions 
of both models are, indeed, equal {6ss — 0.775): not surprisingly, the difference between the two models is in the 
equilibration time-scales, which are considerably faster within the adiabatic approximation. We can therefore conclude 
that the adiabatic approximation is valid on long dynamical timescales. One can also verify analytically that 7c is 
the same for both models: both the adiabatic and non-adiabatic models will undergo a phase transition for the same 
values of the relevant parameters, kc, £, and kg- (The phase diagram for the full N = 2 model is shown in figure |6|. 

E. Introducing Spatial Dependence: Chiral Phase Transition 

Spatial inhomogeneities in the concentrations of the many reactants will have a strong impact on the evolution of 
chirality. As we remarked earlier, the analogy with ferromagnets suggests that spatially extended domains of opposite 
chirality will form and compete for dominance (H [10]. In what follows, we will introduce spatial dependence to the 
full N = 2 model, and study its time evolution as the production rate of monomers from the substrate decreases in 
time. 
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where k is the 
10^''^m^s~^ for 



The advantage of imposing the adiabatic approximation to the N = 2 model is most obvious once spatial dependence 
is introduced: the system is elegantly reduced to an effective scalar field theory where the field determines the 
net chirality in a given volume [1], playing the same role as the magnetization in an Ising ferromagnet. However, 
as we demonstrated in the previous subsection (see figure [2]), the adiabatic approximation changes the effective 
equilibration time-scale, thus affecting the overall evolution of the reaction network. We thus move on to introduce 
spatial dependence to the full N ^ 2 model, following the usual procedure in the phenomelogical treatment of phase 
transition dynamics [T31 113 1 by writing the total time derivatives in eqs. 17 as d/dt — > d/dt — kW^ 
diffusion coefficient. Typical values for k are k = 10~^m^s~^ for molecular diffusion in water and k - 
air. 

The equations are made dimensionless by introducing the time and space variables, to = ^ot f^nd xq = xy^ Xo/k. 
Dimensionful values are then recovered for a particular choice of the parameters kg, Qsi and k. Using the nominal 
values of kg « 10~^^cm'^s~^, Qs ~ 10^^cm~'^s~^, and k for water, one obtains Xq = a/2 x fO^^s^^, which yields to 
t ^ 2.3 X lO-^toy and x ^ 8.5 x IQ-^xom. 

To investigate the dynamical evolution of the net chirality as it undergoes a chiral phase transition, we simulated 
the evolution of the net chirality in a reactor pool by allowing 7 to decrease linearly from 70 = 0.3 > 7c to 7 = 0, 
writing j{t) = 7o(l — t/tf), where tf is chosen so that j{tf) = 0. For the full N = 2 model, the net chirality is given 

by, 



A1 + A2 
S1 + S2 



(19) 



where Ai , A2 , Si , and ^2 were defined in section 111 C 



To generate small initial spatial fluctuations in the concentrations, we coupled the system to an external environment 
via a stochastic spatiotemporal Langevin equation. Rewriting eqs. [iTlwe have. 



-1 (di^ 



° > dt 

-1 fdSi 



A, 



A, 



" > dt 
-1 f dAi 
\ dt 

-1 / dS2 



" > dt 
1 f dA2 
\ dt 



- kW^Si 



k\/^S, 



kSJ^A 



= 1- ^'«(7L + 7fl)'0-Ki/'52 + C(x,t), 

- K (\{iL + 7fl) + 82) v- - Si + e(x, i). 



\{Sf + Al)-S2Si + a^,t), 



-S2A2-A2Si+^{x,t), 



(20) 



where £,{x,t) is a stochastic force with zero mean ((^) — 0) and two-point correlation function (^(x', i')^(x, t)) = 
a^(5(x' — x.)S{t' — t), and measures the strength of the external influence in units of (lcngth)''timc, where d is the 
number of spatial dimensions. 



The equations dictating the evolution of the reaction network (eqs. 20 1 were solved with a finite-difference method 
in a 1024^ grid with St = 0.005 and Sx — 0.2, and periodic boundary conditions. In 2d, this corresponds to simulating 
a shallow pool of linear dimension I sa 200cm. The reactor pool was initially set with a near-racemic distribution 
of monomers and dimers (as would occur in the case that the system started from a small amount of subtrate and 
was permitted to evolve for a period of time) with {9{t = 0)) ~ 10~^. (The angled brackets denote spatial averaging. 



3 In addition to 7 linearly decreasing 
uence such that a^(i — 0) = — 0.02 



(...) = (l/V) / ■ . ■ dV.) The functional dependence for j{t) is shown in figure 
with time, we also imposed a linear time dependence on the environmental inf 
and a{tf = 600) = 0. In other words, 7(a) = 7o(a/ao), thus establishing a relationship between the substrate abihty 
to produce monomers and the external environment. One may think of a{t) as an effective temperature: as a(t) 0, 
the autogenic production of monomers from the substrate ceases. Note that we chose ag to be sufficiently small 
such that the environmental influences do not affect the dynamical evolution of the net chirality. As demonstrated 
by Gleiser and Thorarinson [S], for a above a critical value Oc, chiral symmetry is restored. In the N ~ 2 model 
with adiabatic approximations, = 1.15 in 2d, and = 0.65 in 3d. A recent study of the critical influence of the 
environment on homochirality can be found in ref. |11| . 

The temporal evolution of the spatially-averaged net chirality {0{t)) is shown in figure |4] and two-dimensional 
snapshots of the reactor pool are shown in figure |5] In looking at the evolution of {0{t)), we see that there is a long 
period of time where the net chirality is nearly racemic, even after 7 drops below jc- However, once 7 becomes as 
small as 0.15, larger values of opposing chiralities become possible and we see that the net chirality evolves toward 
the chirally pure phase {0{t)) = —1. 
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t(A„-') 

FIG. 3: Evolution of 7(t) — 0.3(1 — t/600) for the full N — 2 model. The Roman numerals in the figure correspond to the 2d 
snapshots shown in figure |5] 
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t(A„-i) 

FIG. 4: Evolution of the spatially-averaged net chirality {9{t)} for the full N = 2 model with 7 = 7(t) as specified in figure [s] 
The Roman numerals in the figure correspond to the snapshots shown in figure [5] 

In snapshot I of figure [5I 7 = 0.3 and {0{t)) w 10-^ With 7 > 7c the tendency is for the system to remain a 
near-uniform racemate. In snapshots II and III , 7 = 7c = 0.25 and 7 = 0.2, respectively. Even though 7 < 7c, the 
phase space available to 9 is still very small, so we do not yet observe any chiral domains forming. In snapshot IV, 
7 — 0.15 and we begin to observe chiral domains as 9{t, x, y) is permitted to take on larger values. As 7 is decreased 
still further, in snapshots V through IX, the domains coarsen and the magnitudes of both the red {L) and blue (R) 
phases increase. Once domains of both chiralities are present, the curvature pressure of the walls will begin to affect 
the dynamics [10 . Here the red (L) phase does not percolate the lattice so it will continue to shrink in size; eventually, 
only the blue (i?) phase will remain and homochirality will be achieved, completing the phase transition. 
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vii vm IX 

FIG. 5: Evolution of net chirality for the full N — 2 model, with 7 = 7(a), where a is a measure of environmental influences. 
Each snapshot depicts the net chirality 6{x, y) at a time specified in figure |4] Note that 7c = 0.25 is reached in II: for 7 < 7c 
domains of opposite chirality begin to form and compete for dominance. 

IV. THE iV ^ 00 MODEL 

To study the full set of polymerization equations as they might apply to a realistic biochemistry, one must consider 
the reaction network in the limit of large N . Fortunately, as pointed out by WC, the set of infinite rate equations 
may be reduced exactly to a closed system of only a handful of equations, significantly reducing the complexity of the 
problem. 
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To investigate the infinite set of rate equations for our model, we introduce the new variables 

n=2 n=2 

CXD OO 

[U] = J2{LnRl], [T] = J2{RnLi], 
n=2 n=2 

oo oo 

[G] = 5^n[L„], [i?]^5^n[i?„]. (21) 



n=2 n=2 



The reaction network may be written succinctly in terms of these new variables. Setting kj/ks = 1, the entire network 
reduces to (for oo) 



d\S] 
dt 



Qs - {SL + eR)[S] - kc[S]{[L] + [R]), 



^[^^^ - {eL + kcmS]-2ks[L,]([L,] + [R,] + [R] + [L] + l{[U] + [T]) 



'^^^'^ - {eR + kc[R])[S]-2ks[Ri]([L^] + [R^] + [R] + [L] + h[U] + [T]) 



dt 
d[R] 



m 

dt 

dm 

dt 
d[G] 

dt 
d[D] 



ks[L^]^-2ks[Ri][L], 
ks[Ri? - 2ks[LM 
2ks[Ri][L] - ks[Ri][U], 
2ks[L,][R] - ks[Li][T], 
2ks[L,]^ + 2ks[L,][L] - 2ks[Ri][a], 

2ks[Ri]^ + 2ks[Ri][R] - 2ks[L,][D]. (22) 



dt 

Similar to the procedure for the reduced model, we introduce symmetric and asymmetric dimensionless variables 

<5 = i^)H{Li] + [Ri]), A = {'^)H[L^] - [R,]), 
Qs Qs 

^ = (^)^(m + [R]), V = {^^)H[L] - m, 

Qs Qs 

Qs Qs 

■p = i^^)H[G] + m, p = (^)i([G] - m, (23) 

Qs Qs 
along with the dimensionless variable ip = {^f)^[S] describing the substrate. The network in terms of these new 
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variables becomes 



_ 1 d-tjj 1 ^ 







f -fR)Tp + M - s{s - 


f AA + 


idA 
d^ 




- 1r)'>P + KipV - A{S -f 


-M + 


, dA/ 


1 


A')-^-{SAf-Av), 




Tt 


1 


^{Srj-AM), 




^dM 




-Av)-^iSM-Afi), 


, d/i 
^0 d^ 




AAf) - ^{S^L-AM) 




..idT' 
^0 dF 




A^) + ^{SAf + Av) - 


lisv 


^0 dt 




{Sv + AM)-\{Sp- 


AV). 



(24) 

The parameters 7^, 7_r, and k are as defined in section [ill C| 

Note that in terms of these new variables the net chirality is simply 

As shown in figure [6] it is possible to construct a phase diagram for the steady state value of the net chiral asymmetry, 
6ss, by obtaining the equilibrium solutions to eqs. [24] for varying 7. We note that in the N —> 00 limit the only 
possible steady-state solution for 7 = is racemic (not shown in figure [6]): there are, however, runaway solutions 
that, in spite of never reaching a steady state, still reach homochirality. This can be seen from eq. [25| even though r] 
and Af keep growing (the runaway behavior), their ratio approaches unity as A and S approach small values. As WC 
noted, in this case the average polymer length goes to infinity, not a very desirable feature. This supports including 
autogenic reactions in more realistic models. But, as we shall see, these reactions cannot be too efficient or they will 
lead invariably to a racemate, also not a desirable feature. 

In figure [6| we show the phase diagram for the N — 2 and the N 00 models. We have shown above that 
7c|Ar=2 = 0.25. The critical value 7c |iv— >oo= ^ was found empirically: it is the only value of 7 that yields a set of 
solutions which are all racemic and contains no imaginary solutions. (For 7 < 7c chiral solutions exist and, for 7 > 7c, 
there exist imaginary steady-state solutions. Only at 7 = 7c we should see no chiral and no imaginary steady-state 
solutions) . 

The average polymer length is defined as. 



where V and M where introduced in eqs. 23 As mentioned previously, for 7 = the average polymer length diverges. 
However, for 7 non-zero. Lav takes on finite values. The steady-state average polymer lengths for 0.01 < 7 < 0.7 are 
shown in figure |7] For racemic solutions. Lav = 3. This is consistent with the results of figure [7] as 7 approaches 7c, 
the system approaches a racemate and the average polymer is a trimer. 



V. SPATIOTEMPORAL DYNAMICS OF POLYMERIZATION: MAXIMUM POLYMER LENGTH AND 

CHIRAL ASYMMETRY 



In the previous sections, we have alluded to a relationship between the maximum polymer length in the reactor 
pool, N, and the steady-state value of the chiral asymmetry, 9ss- In particular, we ask the question: how does the 
presence of a monomer-producing substrate with rate 7 affect the reactor pool's ability to achieve homochirality for 
different models with maximum polymer length A^? In other words, for a given 7, is it easier or harder to achieve 
homochirality with increasing polymer length? 
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FIG. 6: Phase diagram for the steady state chirahty 0ss as a function of the parameter 7 in the hmit TV — > 00 (black squares), 
where 7 = 71, = jr. Note that jc |jv^oo= \/2/2. Also included is the phase diagram for N = 2 (grey squares), with 
7 = 71, = 7_R, for comparison. The continuous line is the fit 6ss = (1 — ilc/'y)'^)^^'^ 
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FIG. 7: Average polymer lengths as a function of 7 for A'' 00. SS indicates that the values are for the steady state solutions. 



The phase diagrams for N = 2 and A'^ ^ 00 of figure |6] show that, in principle, systems with large N achieve 
homochirality for a wider range of parameter space. Inspecting figure [6) it can be seen that for A^ ^ cx) significant 
homochirality, 9 > 0.9, is achievable for relatively large values of 7, 7 < 0.4, as compared to 7 < 0.06 for N — 2. 

To investigate the spatiotemporal dynamics of the reaction network for different N, we used a finite-difference 
method on a 1024^ grid with St = 0.005 and Sx — 0.2 and using 7 = 0.1 and k = 0.5. For N = 2 with adiabatic 
approximations this entails solving a set of two coupled partial differential equations (cf. 
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11 with d/dt 



eq. 

where we do not evolve 



d/dt — kW^); for the full N = 2 there are five (eq. |17| ; for A^ — > cx) there are seven (eq. 
the V and p equations since they have no bearing on the net chirality and do not affect the dynamics of the other 
7 equations); and for N = 5 there are seventeen coupled partial differential equations (not shown). Each system 
was initially prepared in a near racemic phase with a 10~^ asymmetry in initial monomer concentrations. All other 
polymer concentrations were initially set to zero. 
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FIG. 8: Temporal evolution of spatially-averaged net chirality for N = 2 (grey line), N = 2 with adiabatic approximation 
(dashed black line), N — 5 (dashed grey line), and N ^ <x (continuous black line) for k = 0.5 and 7 — 0.1. (Note that for 
N = 2 with adiabatic approximations, {S{t)) is plotted instead of {6{t)).) 



Figure [S] shows the temporal evolution of the spatially-averaged net chirality, {0{t)), for N — 2 (with and without 
adiabatic approximations), = 5, and N 00. One can see that for all values of N, {6{t)) approaches a constant 
for large times, indicating that steady-state conditions have been achieved. It is also clear that the asymptotic values 
of {0{t)), 9ss, increase with increasing N: for N = 2 (full model), Ogs 0.83; for N = 5, Oss 0.97; and for TV 00, 

To study the spatiotemporal steady-state solutions in more detail, we examined the behavior of systems initially 
prepared in a homochiral phase {{0{t = 0)) = 1) for different values of 7 for models with several values of maximum 
polymer length N, and for N 00. (Again on a 1024^ grid with St = 0.005 and 5x = 0.2.) That is, we examined 
whether homochirality is a stable equilibrium solution of the reaction network, and how this stability depends on 
both the maximum polymer length and autogenic reaction rate 7. The results are shown in figures [9] and [TO] From 
figure [9] we can immediately see that homochirality is only strictly stable as — > 00: for all other N, Oss drifts 
away from unity when perturbed. However, even for models with N as low as = 4, Oss ^ 0.8, or within only 



20% of homochirality. For A^ = 2, 7 = 0.3 is larger than the critical value (7c |Ar=2= 0.25, see section III I so the 



system is stable as a racemate. We have verified that this trend is true for any 7 > so long as 7 < 7c |Af^oo= 
Independently of A^, for 7 = all systems achieve homochirality, while for 7 > 7c Iat— >oo racemates are the only stable 
solution. (Recall, however, that in the Noo limit, the homochiral solutions are unstable and unphysical, as discussed 
in section 4.) Therefore, there exists a range 0/7, < 7 < where the net chiral excess is dependent on N. Using 



the nominal values for reaction coefficients introduced in section III C 



Figure [T0| shows that the dependence of Oss on both A^ and 7 satisfies the approximate fit. 



Oss = (1 - (7c/7)')'/' tanh(aA^ - 67), (27) 

where a and b are 7-dependent fitting constants, and we used 7c — 7c |Af^oo= For 7 < 0.5, covering most of the 
range of interest, a = 0.5 and 6 = 3 give quite accurate results. Figure [TO] shows the curves for 7 = 0.3 and 7 = 0.5. 

A related question is whether the increase in chiral excess with increasing A^ is dominated by longer or shorter 
homochiral chains. The answer to this question will shed light on the actual mechanisms that allow for greater 
chiral purity for large A^ systems. In the context of their spatially-independent model without autogenic monomer 
production, WC found that the total chiral purity of all polymer chains is greater than that of monomers. Here, 
we study in detail how rjn changes with increasing n and varying 7. We obtained the steady state values of the net 
chirality of individual polymer lengths for the set of ordinary differential equations governing the evolution of the 
A^ = 5 system. The results are shown in the stem diagram of figure [TT] One can see that for fixed 7, the longest 
polymers achieve the largest chiral excess. On the other hand, it is also clear that homochiral growth is inhibited by 
large autogenic reaction rates. 
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FIG. 9: Spatially-averaged, steady-state value of the net chirality, 9ss, for various maximum polymer lengths A'^, with k — 0.5 
and 7 — 0.3. Since 7 > 7c |jv=2= 0.25, 9ss for = 2 is zero. Note that as A'' increases so does the net chirality, with 9^1 for 
N ^ 00. 
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FIG. 10: Fit of steady-state chirality, 9s3, for models with various maximum polymer lengths A'', with k — 0.5 and 7 — 0.3 
(black curve and circles) and 7 — 0.5 (grey curve and circles). The fits (continuous curves) are given in eq. 
7 > 7c |jv=2= 0.25, 9ss for A'^ = 2 is zero. The same for A'^ = 3 when 7 = 0.5. 



27 Since 



Sandar's concluded that as N increases bifurcation happens more readily due to the increased number of opportu- 
nities for the minority enantiomer to be removed by cross-inhibition |21j . This conclusion was used to explain why 
large N systems have less stringent requirements on the feedback fidelity in order for bifurcation to occur. Although 
we have not discussed the effects of varying fidelity here, presumably this can also explain why, with the racemizing 
pressure created by 7, it is only in the limit N 00 that 6* — > 1. In order for the net chirality to increase, the rate of 
chiral amplification through enantiomeric cross-inhibition and autocatalysis must exceed the losses from racemizing 
pressures. In previous sections, we have described a competition between the racemizing pressure due to 7 and the 
pressures favoring chirality due to joint enantiomeric feedback and cross-inhibition. The ensuing dynamics are not 
only dependent on the relative reaction rates, e, kg and kc, but also on the length of the longest polymers formed 
through Cl and Cu. The presence of these terms allows for polymers to accelerate the formation of monomers of the 
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FIG. 11: Steady-state values for the chiral asymmetry ?7n(t) for N = 5, with 1 < n < 5. The stationary value of 9 is also shown 
at the far right. Larger polymers achieve a higher degree of chiral purity. Also, the net chirality per polymer is suppressed by 
increasing 7. Note that for t;„ = 1 some points overlay with others. 



same chirality while inhibiting the formation of polymers of the opposite chirality. Due to the increased opportunities 
for removal of the minority enantiomer and for enantiomeric autocatalysis provided by large N systems, it is therefore 
only in the limit of large N that the enantiomeric pressures can overcome the racemizing pressure of a (relatively) 
large 7. 



VI. CONCLUSION 



We know little of the prebiotic conditions that led to first life in early Earth and even less about other possible 
life-bearing planetary platforms in this and other stellar systems. The first amino acids and sugars may have been 
formed here, or may have been fed from outer-space; in either case, there may have been a net initial chirality in 
their concentrations, or the initial conditions may have been racemic. The substrate (or substrates) may or may not 
have supported autogenic production of monomers, with or without a chiral bias. Given these uncertainties, it is 
important to analyze in detail the polymerization dynamics of different reactor pools in order to compute the final 
net chirality produced. Our impetus for introducing autogenic (non-catalytic) monomer prodution to the model has 
been to simplify the initial conditions, thereby allowing for a truly bottom-up approach to the origin of prebiotic 
homochirality. We have found that this mechanism provides for a very rich model of chiral evolution whereby the net 
chirality is dependent upon the autogenic reaction rate parameters 71, and through the length N of the longest 
polymers formed. We have shown that the presence of these terms allows for a chiral-symmetry breaking phase 
transition to occur: for symmetric autogenic production, where 7l = 7_r = 7, we found that, for 7 > 7c (-/V), where 
7c(-/V) is an A^-dependent critical value, racemic solutions are the only equilibrium solutions. In other words, efficient 
autogenic production of monomers from the substrate creates a racemizing pressure that overwhelms the tendency 
toward homochirality from autocatalysis with enantiomeric cross-inhibition: if this production is too efficient, the 
average polymer length is strongly suppres sed (c/. figure[7]), and homochirality cannot be achieved. For 2 < N < 00, 



we found 0.25 < 7c < v2/2. Using eq. 16 we can express this bound in terms of the various reaction rates in 
the model. Adopting a similar procedure, equivalent bounds can be derived for any set of polymerization reactions 
featuring autogenic production terms. 

Our results imply that although in the presence of moderate autogenic monomer production homochirality is 
achievable starting simply from a substrate, it is preferably achieved for reactor pools allowing the formation of large 
polymers. We found that this result can be expressed by an approximate relation for the steady-state value of the 
net chirality, 6ss ~ tanh[7V]. This has implications for studying autocatalytic polymerization networks that display 
chiral symmetry breaking in the laboratory: if non-catalytic production of monomers is active, only systems which 
allow for the formation of large N will support bifurcation toward significant chiral purity. In principle at least, it 
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should be possible to test the onset of chiral symmetry breaking in autocatalytic systems such as the Soai reaction 
by varying external parameters affecting the reaction rates. 

Extrapolating from what we have learned thus far in our solar system and elsewhere, the staggering diversity of 
potential life-bearing planetary platforms spread throughout the galaxy implies in an equivalently staggering diversity 
of prebiotic conditions that possibly led to chiral life. Searching for general, model-independent results and trends is 
thus of the utmost importance. 
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